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Abstract 

Prompted by recent results on Susy-U(N)-invariant quantum mechanics in the large N limit by 
Veneziano and Wosiek, we have examined the planar spectrum in the full Hilbert space of 
U(N)-invariant states built on the Fock vacuum by applying any U(N)-invariant combinations 
of creation-operators. We present results about 1) the supersymmetric model in the bosonic 
sector, 2) the standard quartic Hamiltonian. This latter is useful to check our techniques 
against the exact result of Brezin et al. The SuSy case is where Fock space methods prove to 
be the most efficient: it turns out that the problem is separable and the exact planar spectrum 
can be expressed in terms of the single-trace spectrum. In the case of the anharmonic 
oscillator, on the other hand, the Fock space analysis is quite cumbersome due to the presence 
of large off-diagonal O(N) terms coupling subspaces with different number of traces; these 
terms should be absorbed before taking the planar limit and recovering the known planar 
spectrum. We give analytical and numerical evidence that good qualitative information on the 
spectrum can be obtained this way. 
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1. Introduction 

Interest in the planar (topological) expansion in QFT has always been alive since the work of 
't Hooft |tH74] and Veneziano |Ven74t IVen76| (see also [CMV75j ) . A recent revival of interest 
has been triggered by a paper by Veneziano and Wosiek [VW0_6b] (hereafter VW) followed by 
[VW06cl IVW0 6d] . about a class of super symmetric quantum mechanical models, the simplest 
among them being defined by 

(1) H = {Q,Q^}, Q = j r (p(a + ga 2 )) 

where / and a are NxN operator valued Hermitian matrices, with standard canonical (anti)com- 
mutation relations; the model reveals a rich structure in the spectrum, in particular a peculiar 
duality b — > 1/6, where b 2 = g 2 N is 't Hooft's parameter. The approach to the planar limit is 
based on the Fock space representation of the Hamiltonian, which is for the time being the only 
technique available in cases where a change of variable to [/(iV)-invariant "radial" coordinates 
is not applicable [BIPZ78, MO80J. The results presented by VW are restricted to a peculiar 
subspace of the full sector of U(N) invariant states, namely the subspace TiS Q ^ spanned by the 
vectors which are obtained by applying a "single trace operator" to the Fock vacuum 

(2) Tr(a tn )|0), n = 0,1,2,... 

This subspace is considered "dominant", in the planar limit, with respect to other states built 
over the vacuum by applying any product of invariant operators. The Fock space approach to 
the large N limit was considered in [Tho79| . where the issue of which kind of operators would 
leave the subspace of single-trace states invariant in the planar limit was addressed. The VW 
Hamiltonian ([I]) enjoys this property, while the anharmonic oscillator Hamiltonian 

(3) H= I T r(p 2 + g 2 ) + ^Tr(g 4 ) 

does not. It follows that in this latter case it is necessary to analyze the Hamiltonian in the full 
Hilbert space of U(N)— invariant states. In the SuSy case it is clear that restricting to TC^ is 
legitimate in the planar limit, but one misses in this way a big portion of the spectrum, and, 
moreover, to estimate the 1/N corrections one has to study the full Hilbert space. 

The multiplicity of states makes the calculation of the spectrum much harder. To be specific, 
at a given number of quanta, i.e. considering a product of n creation operators on the vacuum, 
there exists one state in 7i^> but there exist p(n) states of the kind 

(4) Tr(a tAl ) Tr(a fA2 ) . . . Tr(a fAm ) 1 0) , m<n 

where A is any partition of n and the function p{n) counts them all. It is known since Hardy 
and Ramanujan that p(n) grows exponentially (p(n) = 0(exp{K^n}), K = 11^2/3 ), and this 
is the main obstacle in a purely numerical approach. For instance, while it is easy to diagonalize 
the VW Hamiltonian corresponding to a maximum number of quanta of thousands or even more, 
the mere fact that p(200) «4x 10 12 (as it was first computed exactly by MacMahon |Har02j) 
puts a severe limit to the level one can explore numerically. 

Our main result (Sec. 2) will be the following: for the SuSy VW Hamiltonian in the bosonic 
sector the Hilbert space of U(N) invariant states splits into an infinite number of subspaces 
H ( - m \m = 1,2,3,..., which are left invariant by the Hamiltonian in the planar limit. These sub- 
spaces are characterized by states with a given number m of trace factors in Eq.([3]); so to speak 
the number of traces is a good planar quantum number. As a consequence, the Hamiltonian is 
recognized to be separable into the sum of commuting operators, all unitarily equivalent to the 
single trace operator studied by VW. Essentially the Hamiltonian is analogous to an isotropic 
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multidimensional harmonic oscillator, where the one-dimensional operator is given by VW's op- 
erator. This fact has several consequences, for instance if the 't Hooft parameter b 2 = g 2 N is 
grater than one, then all levels are infinitely degenerate. 

The corrections of 0(1/N) break this degeneracy since they couple the subspaces ii~( m ) to- 
gether. We do not have exact results on this breaking mechanism, only numerical evidence. It 
appears that while higher states' degeneracy is lifted, the first level (the zero mode) keeps its 
infinite degeneracy even after the 1/N corrections are taken into account. 

For a second class of Hamiltonian operators, which typically do not leave the Fock vacuum 
invariant, as is the case for the anharmonic oscillator solved in [BIPZ78], the analysis is more 
involved. We give evidence on how the exact result, derived by introducing "radial variables" 
for the iV-dimensional matrix x, can be recovered numerically by reabsorbing the off-diagonal 
0{b 2 N) terms which deny the possibility to have a straightforward planar limit. Our conclusion 
is that Fock space methods are not particularly convenient for this kind of models, but they can 
still be useful to explore the spectrum numerically in those cases where a BIPZ technique is not 
applicable. 

2. Susy quantum mechanics in the U(N) invariant sector 
Let us introduce a bit of formalism. Let 

(5) | A) = |Ai, A 2 , . . . , A m ) = M x Ma^)Jr(a^) . . . Jr(a^)\o) . 

The right-hand side is symmetric in the exchange of A's, hence the left-hand side will always be 
reduced to a normal form subject to the condition 

(6) Ai > A 2 > . . . > A m . 

Where confusion may not arise, we shall denote the basis vectors simply by In, m, A) or simply 
| A) where A is an m-partition of n £ Z + , a non-negative integer; we write |A| = n or A h n if A 
is a partition of n as it is customary in the mathematical literature. The "one-trace-states" of 
Ref . [VW06bj are special cases, namely n, 1, {n}). 



If we refer our states to a basis of coherent states \z}, such that 

(7) Qij\z) = Zij\z) 

we realize that the states n, m, A) realize a basis in the space of symmetric functions in 
the eigenvalues of z. Other basis, such as Schur's, can also be used where convenient (see 
e.g. [Led77l ISta99] ): in Sec. 3 we shall turn to Schur's basis which coincides with the free 
eigenbasis (for a detailed description of this formalism, see |Sto90] ). One should take into 
account the fact that our basis is not orthogonal. This is not going to raise any problem, the 
only consequence being that the matrix representing the Hamiltonian is not symmetric, but still 
unitarily equivalent to a Hermitian one. It is however important to normalize the states at least 
to leading order, to avoid having unnatural dependence on N in the matrix elements. So we 
decide to fix the normalization factor to take into account the leading power in ./V: 

(8) M\ = AHAI/2 _ 

(this is just the scaling with N, the exact normalization is known, but not actually necessary, see 
Appendix 2). Now we compute the matrix elements of the Hamiltonian introduced in Eq.([T]), as 
it is reduced to the bosonic sector. For details about the algebra involved the reader is referred 
to |VW06al IVW06bj . We are led then to compute 



(9) (jr{a ] a) + gJr{a) 2 a + a^a 2 )+g 2 Jr{a) 2 a 2 )\ \\) = J^Wa',A |A'> 
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By applying the commutation relations 

(10) Wij,%l\ =S H S jk 

one must bring all annihilation operators to the right eventually reaching the vacuum. This is 
easily done symbolically using a language like Mathematica® , but it is clear that the operation 
is going to exhaust the capabilities of your computer as soon as the number |A| reaches ten or a 
little more, because the multiplicity of terms generated. A way out to this practical limitation 
is provided by the following analytic result (notice that the notation Xj means that the term Xj 
must be deleted from the list): 

Theorem 1. 

i) 

Tr(at 2 a) | Ai, A m > = N 1 / 2 £ Xj \X X , X, + 1, A m ) 

i 

a ) 

Tr(a t2 a 2 ) |Ai, A m ) = ^ Xj y~) |Ai, -, Xj - s, X m , s) 

l<j<m s=0 
+ 2 ^2 Aj A; | Ai , Xj, A;, A m , Xj + A;) 

Hi ) 

T^^lAi,..,!™) = N- l l 2 £ ^ Yl \h,-Aj-s,...,X m ,s-l) 

j s=l 

+ 2 N~ l l 2 Xj Xi | Ax, Aj, A;, A m , Xj + A/ — l) 
Proof. The proof follows essentially the ideas introduced in [Tho79| : 

i) having just one annihilation operator to commute, the result is given by the sum of terms 
obtained by applying the commutator to each factor in |A), and the result coincides with that 
one has for single trace states. The number of quanta on the right increases by 1, hence the 
power of iV which takes into account the normalization. The resulting kets should be put in the 
standard form of Eq.© after the operation has been completed. 

ii) The two annihilation operators can be made to commute with the same factor a^j , provided 
Aj > 1, in which case we have a splitting into two traces or with two different factors, in which 
case we have a merge into one single trace. A factor 2 is due to the fact that the same pair (j, /) 
can be coupled in two different ways to a 2 . The term with s = produces a factor N and this 
is the dominant term in the result. As before the ket arguments on the right-hand-side should 
be rearranged in order to have a non-increasing m-tuple. 

in) As before, the two annihilation operators can be commuted with the same Tr(a^ Aj ) term, in 
which case we have a splitting into two traces; when s = 1 we get a factor which gives the 
dominant contribution. Otherwise the two a's commute with two distinct factors and we have 
a merge into a single trace. □ 

Notice that in case (ii) the operator has a g 2 factor, which means that the dominant term 
O(N) produces a finite result in the planar limit, the other being depressed by a factor 1/Af. 
In the other two cases the factor g combines with the N 1 / 2 factor to give again a finite limit. 
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Other terms provide off-diagonal corrections 0(1/N), which can be taken into account to get 
the planar expansion. 

Finally we have the matrix elements independent from g = 0, that is Tr(cjt a) In, m, A); this is 
a diagonal operator, with eigenvalue n = |A|, independently from the number of traces m. 

To get an idea of the full matrix, consider the portrait in Fig.l (left), which includes all non- 
vanishing matrix elements, including the 0(1/N) terms, for n < 16 (914 states) Much simpler 




Fig. 1. The sparsity pattern of TL up to A h 16 (left); the same on the right after 
taking the planar limit. 

is the planar matrix, the finite part at N = oo, which is depicted in Fig.l on the right. Both 
patterns are replicated with no substantial variations for higher values of n (we checked up to 
n = 22). 

In spite of the apparent complexity of this matrix, which would discourage from pushing to 
higher partition numbers - remember the exponential growth of p(n) - it turns out that the 
planar limit is under full analytical control. This is due to the following observation: the finite 
matrix elements in the planar limit are all generated when the Hamiltonian operator interacts 
with a single factor Tr(a^ Aj ), and the result is built by the sum of all these contributions. This 
means that in the planar limit the number of traces in |A) is unchanged by the application of 
H, in other words the Hilbert space splits into invariant subspaces 7^( m ) characterized by a fixed 
number m of traces and the planar Hamiltonian is given by the sum of commuting operators 
Hj, each of them acting on the j — th index in the partition. Moreover each Hj is identical 
to VW operator in the bosonic sector, that is the Hamiltonian restricted to single-trace states. 
Hence we have 

Theorem 2. The planar spectrum of the VW SuSy operator in the full Hilbert space of U (N) 
invariant states, as defined in Eq. ([9]), is given by the union of the spectra in each subspace 
TC^ m \ m = 1, 2, 3, namely 

(11) E ni ,...,n m = Eff + ... + £«, m > n 2 > ... > n m 

where En^ are the eigenvalues of the Hamiltonian restricted to the single trace states as given 
in [VW06b| . 
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Proof. The result follows from the property that the number of traces is left invariant in the 
planar limit; the constraint on the quantum numbers stems from the symmetry under permu- 
tation of the trace factors. (To make the paper more self contained, we give a derivation of the 
VW spectrum in Appendix 1). □ 

Notice that the block-diagonal structure of the planar matrix can be made manifest if we 
rearrange the partitions in order of increasing number of traces, as shown in Fig. 2. 




Fig. 2. The rearrangement of Ti (left) and its block-diagonal structure in the 
planar limit (right). 

We tested the result by diagonalizing each diagonal block. In this case we can explore larger 
matrices and the result is checked to high accuracy. 

Some comments on the structure of the spectrum are here in order. In the weak coupling limit 
the eigenvalues En ~ (1 + b 2 )n, (n > 1), so the spectrum tends to be highly degenerate, like 
the spectrum of the union of m-dimensional harmonic oscillators for m > 1. Clearly, the states 
with m traces set in only starting with at least m quanta, which implies that the multiplicity of 
levels is not infinite, but it is steadily growing. See Fig. 3 for the case b = 1/4, where the typical 
p(n) multiplicity is manifest. 

The spectrum is dramatically different at strong coupling. As it has been shown in [VW06bJ, 
there is a duality between weak and strong coupling, given by 

(12) + l)/b = b(E n %(l/b) + 1) , (b 2 = g 2 N < l,n > 1) 

Moreover a new zero energy state arises at b > l,E[ 1] (b) = 0. Since the VW spectrum determines 
the whole planar spectrum, the duality relation extends, with slight modifications, to the general 
case. However the structure of the spectrum at strong coupling is deeply changed, namely each 
eigenvalue is infinitely degenerate. This is due to the fact that in the relation (jlip we may 
choose n-i = n% = ... = n m = corresponding to the zero modes and therefore there is an 
eigenvalue EfP in each invariant subspace 7i^ m \ This fact extends to all other eigenvalues: 
taking = = ... = n m = 0, m > 2, one gets an infinite number of eigenstates with the same 
eigenvalue E$ + E^ and so on. 
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This property can only be approximately displayed in a numerical plot, since the diagonaliza- 
tion process implies a truncation on the partition number, however a clear plateaus formation 
can be seen in Fig. 3 corresponding to b > 1 (the plateaus would be infinitely long in the exact 
solution). 




20 40 60 80 20 40 60 



(N = 40) (N=°°) 
Fig. 3. The spectrum of Tt, truncated at A h 22, for N = 40 and N = oo. 



What happens of this infinite degeneracy at finite ./V? We do not have analytic results on 
this matter. However diagonalizing the full Hamiltonian including all subleading terms one may 
hope to get some indications. In Fig. 3 we report the two cases N = 40 and N = oo, truncated 
at A h 22 (a 4507 x 4507 sparse matrix). We may observe that the degeneracy of higher levels 
tends to be lifted by 0(1 /N) corrections, while the ground state appears to be robust against 
the corrections. We are lead to conjecture that the zero modes are present even after switching 
on the subleading terms. 



3. The U(N) symmetric anharmonic oscillator 

The Hamiltonian of Eq. ([3]) has been solved in the planar limit in [BIPZ781 IMO80| . Here we 
want to analyze the problem from the point of view of Fock space. The Hamiltonian deeply 
differs from the SuSy model; for instance the Fock vacuum is not an eigenstate, secondly, and 
more important, the Hilbert subspaces 7i.( m ^ are not invariant. These facts make the Fock space 
approach much more involved than in the SuSy case. Nevertheless we find it instructing how 
the exact result may eventually emerge from the Fock space analysis. 



The Hamiltonian contains 
is studied): 
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three main contributions (see [VWQ6aj where a similar Hamiltonian 



(13) H = H + H 2 + H 4 



H = \N 2 + Tr(ata) + ^Tr( a t 2 a 2 + ... + aV 2 ) 
b 2 

H 2 = —Tr(a t3 a + ... + aV) 
= ^Tr(at 4 + a 4 ) 



where the dots point to all other products generated by (a + a^) 4 , which differ from the terms 
displayed by the order of the operators. It would be nice to have simple recipes like in the 
previous section to compute matrix elements. In this case however the problem is much more 
involved, with 16 quartic operators, some of them containing two annihilation operators or more. 
To get an idea of the work required, just consider matrix elements between states with partition 
number (occupation number) ten, that is a 42 x 42 block on the diagonal: this makes 10584 
matrix elements. If we have to reach a matrix truncation corresponding to a maximum |A| of 
order at least 20, the total number of matrix elements to be computed touches the 10 8 range, 
clearly unfeasible by brute force. But also an approach through a computer algebra system is 
going to face serious difficulties if we want to raise the truncation number above 10, because 
of the proliferation of terms in the reduction of operator products. So we need some kind of 
human-assisted-computer-algebra-approach. 



3.1. Matrix elements by computer algebra at low |A|. We have set up a Mathematica code 
which is able to compute the matrix elements by blindly applying the commutation relations. 
Let us observe that since we initially work with the basis of Eq.(|5]) which is not orthogonal, the 
matrix made of scalar products (A H Lu) must be multiplied on the left by the inverse of the 
metric matrix Aix^ = (A|^) before we can use it to compute the spectrum. This is the true 
representative matrix that we denote by TC 



(14) Hx ll = M^(v\H\ fi 



We start examining the main block on the diagonal of the Hamiltonian in the "partitions basis" 
| A). From the analysis of low dimensional blocks a regular pattern emerges. After subtracting 
a term ^N 2 (l + b 2 ) proportional to the identity operator, we are left with an "almost" diagonal 
matrix, in the sense that the off-diagonal matrix elements are depressed by a factor 1/N. For 
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instance: 
Woo 
Tin 

^Ah3,Ath3 



U2 



b z + 1 + 3f^ 

+ 2(1 + 36 2 ) 

3ft 2 . 

JV 

+ 3(1 + 36 2 ) 

9ft 2 . 

JV 



+ 4(l + 36 2 ) 

126 2 

6/j 2 . 
JV 






U2 



Id 



b 2 + 2(1 + 36 2 



66^ 
TV 

■ 2 + 3(1 + 36 2 ) 

JV 

9& 
JV 

> 2 + 4(1 + 36 2 ) 


JV 







9ft 2 . 

7V 

±& 2 + 3(l + 36 2 )/ 

12b 2 
JV 



+ 4(l + 36 2 ) 

66 2 1 h 2 

JV 4" 





l* 2 





12b 2 

K 

JV 

+ 4(l + 36 2 ) 

JV 







186 2 
JV 

i6 2 + 4(1 + 36 2 )/ 



A part from the obvious occurrence of a multiple of the identity of dimension p(n) with a 
value n(l + 6 2 ) + jb 2 , which obviously generalizes to all higher dimensional blocks, the rest of 
the matrix is not readily generalizable, at first sight. 

The off-diagonal blocks are generated by H2 and H4. Let's examine the first few of them to 
see the pattern (we discard terms 0(l/iV 2 )): 



Ti 



Ah2,^th0 



Ti 



Ah3,/ihl 



Ti\\- 



4, M h2 




^Ah0,/Ltl-2 
HAhl,/xh3 



2b 2 N+% 3b 2 



9b 2 2b 2 N+13% 96 2 



12b 2 

12| 6 2 



18& 46 2 AT+26g 




36 2 



2fj 2 7V+25^ 186 2 



The matrix elements for H4 are similar in that they contain "superleading" terms 0(b 2 N) 
which do not allow for a straightforward planar limit. Unlike the 0{N 2 ) along the diagonal 
which only contributes to the zero-point energy, the big off-diagonal terms cannot be disposed 
of in a simple way and they represent the main obstacle to compute the planar limit. With 
hindsight from [BIPZ78j, we know that the zero point energy has a full expansion in powers of 
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b 2 , which means that contributions coming from the off-diagonal terms are needed to recover 
the exact result. We shall postpone the discussion on this point after we switch to another basis, 
that of Schur's functions. 

3.2. The basis of Schur's functions. It is well-known from the theory of symmetric functions 
(see e.g |Led77[ TSta99j ) that there are many interesting basis in the ring of symmetric functions 
on iV indeterminates, A(iV). Up to now we worked with the "partitions basis" of Eq.(|5]); it's easy 
to relate this basis to the corresponding basis in A(iV). Any vector A can be represented in the 
basis of coherent states where it becomes a symmetric polynomial of the eigenvalues Cl>— >CiV 
of the matrix z 

(15) (z|A> = Tr(^)...Tr(z A -) = ^C^ Etf* 

It is now possible to consider other basis in A to explore the properties of H. One such basis 
is defined by Schur's functions (see |Sto90] for an application). This task can be easily achieved 
by symbolic computer algebra. Let us denote by H the representation of H in Schur's basis. 
The main feature of 7i s turns out to be the following: all blocks along the diagonal are now in 
exact diagonal form; the diagonal matrix elements in the n— th block are now easily identified, 
being given by the already computed n(l + 36 2 ) + \b 2 plus corrections 0(1/ N) which follow a 
regular pattern, e.g. for 1 < n < 6 they are given by 



n 


5E = h l /N x 








1 





2 


(3,-3) 


3 


(9, 0,-9) 


4 


(18, 6, 0,-6,-18) 


5 


(30,15, 6, 0,-6,-15,-30) 


6 


(45,27,15, 9, 9, 0,-9,-9,-15,-27,-45) 



Table 1. The fine structure of anharmonic oscillator to first order in 1/N. 



The pattern may not be clear at first sight, but it emerges immediately if we consider the 
partitions and Ferrers- Young diagrams associated to them, as it is familiar from the theory of 
the symmetric group; namely to each partition Ai > A2 > ... > A m it is useful to associate a 
diagram composed of m lines each containing A& square cells, i.e. 



{5,2,1} 



The formula reproducing the data of Tab. 1 is extremely simple 
(16) 5E X = 3 £ (i-j)b 2 /N. 
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where the pair runs on the cells of the Ferrers- Young diagram of the partition (i is the 
column and j the row index, respectively). For instance at n = 6 we find 



A 
A 



{6}- 
{5,1} 



• 3 x ^ [0 1 2 3 4 5] 
_ ^[01234 



45 



27 



By construction the correction is opposite in sign for conjugate partitions, and it averages to 
zero on every diagonal block, a feature which was already noticed in [MO80J. 

Also off-diagonal blocks greatly simplify in Schur's basis. First of all, since the basis is orthog- 
onal, TL S is symmetric, hence we save half of the computing effort. Secondly, the superleading 
terms turn out to be very regular: all coefficients of b 2 N are either 1 or —1, for instance 

/ 1 1 -1 0\ 
Uf h3 Mh5 = b 2 N I 1 -1 + finite terms 
\ 1 -1 -1/ 

3.3. The role of superleading terms. We now want to discuss how superleading terms of 
0{b 2 N) contribute to the planar expansion and allow in principle to recover the known result. 
At this stage we have an expansion for the matrix representing the Hamiltonian in Schur's basis 
of the following form 

(17) H s = ±A 2 (1 + 6 2 )I + NT 1 +T + N- x T-i + ... 

The term T\ should be taken into account before taking the limit N —* oo. It is clear that 
off-diagonal terms of order N will contribute corrections of order N 2 to the eigenvalues, however 
there should not be any higher order corrections, since there are no such terms in the planar 
expansion; moreover these corrections should be the same for all eigenvalues from what we know 
on the basis of the exact solution. The effect of T\ should then reduce to a mere redefinition of 
the vacuum energy. We can deduce even more: T\ cannot have by itself any discrete eigenvalue 
ti, otherwise, in the strong coupling limit, we would find, by perturbation theory in 1/6, a 
leading behavior b 2 Nt\ which is absent in the exact solution. Hence we are led to argue that 
T\ must be an operator with continuous spectrum which perturbs To by simply shifting the 
zero energy level, but leaving the energy gaps unchanged. An example of this kind comes from 
elementary quantum mechanics. Consider the simple harmonic oscillator perturbed by a term 
—Fx; its perturbative expansion stops at second order for all levels and its only effect is to shift 
all the spectrum by a fixed amount — \F 2 . We argue that this is happening with T\ and this 
can be proven by a direct analysis showing that T\ is unitarily equivalent to a component of the 
position operator in an anisotropic two-dimensional harmonic oscillator (see Appendix 2). As 
in the harmonic oscillator example, it would be natural to change the vacuum state in order to 
re-absorb the term T\. This will not be attempted here and is left for future investigational- 

Let us observe that the general picture we get at this stage, a ground state proportional to 
e(b) N 2 and a splitting of order 1/N given by Eq. (I16p is qualitatively correct, and this result is 
due to the introduction of a convenient basis, provided by Schur's functions. 

3.4. The Hamiltonian matrix from group theory. Even if Schur's basis allows for a better 
starting point in attacking the problem, still the limitations posed by the brute force symbolic 
calculation do not allow us to reach a reasonable truncation of the Hamiltonian. But there is 
another way in which the theory of symmetric functions, or in other words the theory of the 



This idea was formulated by G.Veneziano, private communication. 
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symmetric group, can help in the calculation of matrix elements. Let us denote by \\) s the 
vector corresponding to the Schur's function 



s\{x) 



det(x 



det(x*- , 

while we keep the symbol A) for the states constructed via application of multiple traces 
(Eq.Q). The transformation between the two basis 



(18) 



\ x )s = Y Sx ^) 



is well-known and is directly related to the character tables of S n , the group of permutations 
(see |Led77j . Ch. 4): 

where n = |A| = \pt\, is the dimension of the irreducible representation indexed by \i and x is 
its character. We show now how the knowledge of S can simplify the computation of the matrix 
elements. Let us compute 

Tr(at 4 )|A> 5 = ^5 v Tr(at 4 )|/,> = 

/i M'=MJ{4} 
S ^ S ti'\\ U )s 

(here A U \x denotes the standard set-theoretic union of two lists) . The point is that the character 
table can be easily built to rather large dimensions using the results of Jacobi and contemporaries 
and this allows to break the limit of the brute force calculation. The normalization factor is 
easily found to be 

5<A|/i> 5 = 5 Xfl J] (N + i-j) 
(ij')eA 

where the indices i,j run on the cells of the Ferrers- Young diagram of the partition (i is the 
column and j the row index, respectively), e.g. 





n 


21 


31 


41 




12 3 


{4,2,1}^ 


12 


22 






-1 




13 








-2 



which gives the normalization factor 

5 <{4, 2, 1}|{4, 2, 1}> 5 = N(N + 1)(N + 2)(N + 3)(iV - 1)N{N - 2) 

A similar formula can be found also for the matrix elements of H2 in the Hamiltonian con- 
taining three creation and two annihilation operators. In this way we have built the matrix 
in the (orthonormal) Schur's basis up to A h 21 and presumably this may be improved. The 
matrix can be used to explore the spectrum in the planar limit by expanding in powers of 
1/N and neglecting terms 0(1/N 2 ) or smaller. In particular we can study the crucial question 
about the ground state, which should be compared with the exact result (see Fig. 4 where the 
straight line corresponds to ^iV 2 (l + b 2 )). As a final comment, the introduction of Schur's basis 
appears the most natural choice in view of the structure of the exact solution, a collection of 
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Fig. 4. The ground state energy as a function of b 2 at N = 10. 



N non- interacting fermions (see [BIPZ78J). Schur's functions are also known as Slater- Fock 
determinants, of course. 

Conclusions 

Fock space methods are very efficient for a class of [/(iV)-invariant Hamiltonians similar to 
VW model, where the Fock vacuum is stable. We derived the exact spectrum for the SuSy 
model in the bosonic sector, and this should be easily extended to the one-fermion sector. The 
one-trace sector, while representing a small subspace of the whole Hilbert space, provides the 
building block for understanding the whole spectrum. The b 2 = 1 transition marks a boundary 
between a typical approximate p(n) degeneracy to a regime with infinite degeneracy, which is in 
general broken by 0(1/N) corrections which we have identified. Extension of these techniques 
to non-vanishing fermion number will be considered in a future publication. For the anharmonic 
oscillator, as well as for all other Hamiltonians covered by the exact BIPZ solution in the planar 
limit, Fock space methods are not convenient. Still, by a combined effort of analytical and 
computer algebra techniques, they may prove to be reliable in those cases where an approach 
based on "polar coordinates" a la BIPZ is not available. Our results show that superleading 
matrix elements 0(b 2 N) do not spoil the planar expansion, and the spectrum can be at least 
computed numerically Q 
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Appendix 1 



The spectrum of the SuSy model in the single trace subspace has been computed exactly in 
[VW06b]. We give here an equivalent derivation, essentially algebraic in character. 
The matrix to be diagonalized is the following 



(n\H\m) = (l + b 2 (l-5 nl )) 

n &nm + b \J m(m + l)5 n ,m+i + b \J n{n + l)<5 mjn +i . 
Another matrix, differing only at the (1,1) element turns out to be easily diagonalizable: 



(n\Tio\m) = (1 + b 2 )nd nm + by / m(m + l)5„, m +i + b v / n(n + l)<5 m , n+ i . 
To see this, we first change to another representation 

its action being given by 

Tt (j) n = (1 + b 2 ) ncj) n - b 2 (n + l)0 n +i - (n - l)</> n -i 

Lemma 3.1. Ho leaves all the subspaces = {4> n = Nn k + 0{n k ~ 1 )}, (k = 0,1,2,...), 
invariant. 

Proof. By inspection, the terms n k+l exactly cancel and the n k term is multiplied by (1 — b 2 )(k + 
1). ' □ 

Theorem 3. H.q has discrete spectrum given by E n = |1 — b 2 \n, (n = 1,2, 3,...). 

Proof. In each finite-dimensional subspace the matrix is upper-triangular, hence its eigen- 
values can be read off the diagonal. The eigenvectors are given by the basis of orthogonal 
polynomials w.r.t. the discrete measure 



oo 

|2 



^2nb 2n \0 n 



n=l 



which are known as Meixner polynomials |Sze75j . They can be defined by 



( 1 - 6 )V^) 2i?l( " A; ' n+L2 - 1 

and satisfy the orthonormality property 

oo 

J2nb 2n 4> { n k U { n h) =S K 

71=1 



□ 



The rank-one-perturbation property [GL96] gives now the spectrum of M in the F = sector. 
The eigenvalue equation reads 



^(k + l)(l-b 2 )-E 



k=0 

where v% is the first basis vector in the representation of ([1]), and we find 



~ jfe6 2(fc-l)( 1 _ & 2)2 

b > — = 1 

t[ (l-V)k-E 



14 



hence, by setting E = (1 — b 2 )e, one gets 



e - A; 1 - b 2 

k=l 

The equation can be re-expressed in terms of hypergeometric functions: 

A(e) = l + -i^- 2 Fi(l,l -e,2-e,b 2 ) 
e — 1 

which is equivalent to what has been derived in [VWOGbJ by another method. In our derivation 
the structure of A(e) arises naturally via the "rank-one-perturbation" theorem. 



Appendix 2 

We want to discuss the structure of the superleading term in Eq. (|17p . Since we work at N = oo 
we can study the representation of T\ in the partitions basis, which in this limit is orthogonal. 
The leading behavior of the normalization factor can be easily expressed by representing a 
partition in terms of "composition", i.e. 

A = {l ri 2 r \...k Tk ] 

the integers Tj denoting the multiplicity of j in the partition of n = Y^j j r j ■ We find 

f^ 2 = (\\\) * NWHp rj\ 
j 

Now let's apply Tr(a^ 4 ) to the normalized AY We find 

Tr(at 4 )| A) = Tr( a t 4 ) \r u r 2 , . . . , r k ) « ^|n, r 2 , .., r 4 + 1, ...) = 2N 2 Vu+T | A') 

where A' = Au{4}. The operator Tr(a^ 4 ) when multiplied by jb 2 /N gives rise to one component 
of T\ which couples states in the form 

A -> A U {4} -> A U {4 2 } -»■ A U {4"} . . . 

with matrix elements which are identical to those of a single creation operator. The action of 
Tr(a 4 ) is simply obtained by transposition; if a partition does not contain 2 or 4 it is annihilated 
by Tr(a 4 ), which means, of course, that the resulting matrix elements are not O(N), but finite. 
The analysis of the remaining operators in T\, coming from H2, is much more involved, but 
the structure can be identified by using computer algebra. Eventually we get the following 
representation 

T 1 |r 1 ,r 2 ,...,r 4 ,...> = V2b 2 N (y/r^\r h r 2 - 1, u, ...) + \/(r 2 + 1) |n, r 2 + 1, r 4 , ...» 

+ \b 2 N (y/rl \ri,r 2 , ...,r 4 - 1, ...) + \/r 4 + 1 \n,r 2 , ...,r 4 + 1, ...)) 

Notice that in this way we can identify an infinite sequence of subspaces invariant under the 
action of T\: each subspace is labeled by a partition characterized by r 2 = r 4 = 0. This 
is a "vacuum" state for T\. The action of T\ builds a subspace isomorphic to the Hilbert 
space of a two-dimensional harmonic oscillator and the action of T\ in all subspaces is always 
the same. The harmonic oscillator frequencies can be read off the diagonal part To, namely 
(tt»i, uj 2 ) = (2, 4)(1 + b 2 ). Hence we find that indeed the only effect of T\ at this level consists in 
a shift of the whole spectrum proportional to b 2 N exactly as it happens to a simple harmonic 
oscillator under perturbation by a term oc (a + cr). 
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